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Abstract 



We discuss the structure of infrared and ultraviolet singularities in on-shell 
QCD and supersymmetric QCD amplitudes at one- loop order. Previous 
results, valid for massless partons, are extended to the case of massive partons. 
Using dimensional regularization, we present a general factorization formula 
that controls both the singular e-poles and the logarithmic contributions that 
become singular for vanishing masses. We introduce generalized Altarelli-Parisi 
splitting functions and discuss their relations with the singular terms in the 
amplitudes. The dependence on the regularization scheme is also considered. 
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1 Introduction 



Hard-scattering processes that lead to multiparton final states with heavy particles are 
important for present physics studies within and beyond the Standard Model. Their im- 
portance will further increase at future high-energy colliders (see, e.g., Refs. [|l], |2| and 
references therein). Reliable theoretical predictions for these processes require the eval- 
uation of (at least) next-to-leading order QCD corrections, independently of the nature 
(QCD, electroweak, SUSY) of the interaction that produces these processes at the lowest 
perturbative order. The computation of QCD radiative corrections to cross sections that 
involve heavy partons are certainly more complicated than in the case of massless partons, 
as known since the first NLO calculations of heavy- quark production in hadron collisions ||. 

Most of the available techniques (see, e.g., the list of references in Sect. 4 of the QCD 
Chapter of Ref. 0) to perform NLO calculations require the analytic evaluation of the 
infrared (soft and collinear) and ultraviolet singularities of the one-loop amplitudes for the 
corresponding process. The knowledge of the singularity structure in a general (process- 
independent) form is thus useful for several reasons. Loop calculations are always very 
cumbersome and error-prone, therefore, knowing any general property of the result is a 
good check on the calculations. The knowledge of the singularities of the loop amplitude 
can also be used as a possible tool to split the loop calculation into a divergent part, which 
is known in analytic form, and a finite part, to be evaluated by numerical methods. In the 
case of infrared-safe observables the infrared singularities of the real and loop corrections 
must cancel, therefore, knowing the singularity structure without performing the explicit 
loop computation helps devising process-independent calculational schemes at NLO. 

The singular behaviour of QCD amplitudes involving massless partons is completely 
known at one-loop order and the coefficients of the divergent l/e n poles for n =2, 

3 and 4 are also known for two- loop amplitudes in the most general case 0. Despite the 
numerous one-loop calculations involving massive particles that have been performed so far, 
to our knowledge the general form of the singularities in amplitudes with massive partons 
has not yet been presented even at one-loop order. 

In this paper we discuss the structure of the singularities of one-loop QCD (or QED) 
corrections to on-shell scattering amplitudes that involve both massless and massive par- 
tons. By singularities we mean not only the terms that diverge in absence of regularization, 
but also those that become divergent in the limit of vanishing parton masses. With in- 
creasing collision energies the ratios of the particle masses to the kinematic (Mandelstam) 
invariants tend to zero. Thus the amplitude contains logarithms of such ratios that be- 
come large and whose explicit control is important for numerical stability. This feature has 
particular relevance for electroweak physics and one-loop QED corrections where collinear 
singularities are physically regularized by the small values of the light-fermion masses. 

We present a universal (process-independent) factorization formula that embodies the 
structure and the coefficients of all these singular terms in the one-loop amplitudes. We 
use dimensional regularization in d = 4 — 2e dimensions for both ultraviolet and infrared 
divergences, and we thus give the coefficients^ of the corresponding poles 1/e 2 and 1/e. The 

* These coefficients can also be obtained by using the general results of Ref. S. 
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factorization formula also explicitly exhibits all logarithmically-enhanced terms In 2 m and 
In m that become singular for vanishing mass m of the massive partons. Furthermore, we 
can have full control of all the terms that are not analytic in the massless limit, because the 
factorization formula also contains the constant (when m — > 0) terms that originate from 
the non-commutativity of the limits m — > and e — > 0. Therefore, in the limit m — > 0, the 
finite (when e — > 0) one- loop contribution that is not included in the factorization formula 
tends smoothly to the finite contribution of the corresponding massless amplitude. 

There are two ways to derive the singular terms. They can be obtained either by ex- 
plicitly performing loop calculations or by using unitarity. We use the latter method. It 
amounts to exploiting the fact that the infrared singularities of the real and loop correc- 
tions cancel each other in the computation of infrared-safe observables. Thus, the singular 
structure of the loop corrections can be derived from that of the real corrections. Our 
process-independent calculation of the real corrections is based on the dipole subtraction 
formalism ||. In particular, we have explicitly extended the results of Refs. ||[| to the 
general case of partons with arbitrary masses. In this paper we present only the final re- 
sults on the singular behavior of one-loop amplitudes. Complete details of the extended 
formalism will be presented elsewhere. 

After fixing our notation and conventions in Sect. we present our results on the 
singular terms of QCD and SUSY QCD amplitudes in Sect. |3|. The singularities produced 
by QED corrections to electroweak and QCD amplitudes can be obtained from the QCD 
results as a special case. In Sect. f| we relate the singularities to individual terms in 
generalized Altarelli-Parisi functions. Section [| contains our summary, and the Appendix 
provides a list of generalized Altarelli-Parisi functions. 



2 Notation 

2.1 Dimensional regularization 

In the evaluation of loop amplitudes one encounters ultraviolet and infrared divergences that 
have to be properly regularized. The most efficient method to simultaneously regularize 
both types of singularities in gauge theories is to use dimensional regularization |10| . The 
key ingredient of dimensional regularization is the analytic continuation of loop momenta 
to d = 4 — 2e space-time dimensions. Having done this, one is left with some freedom 
regarding the dimensionality of the momenta of the external particles as well as the number 
of polarizations of both external and internal particles. This leads to different regularization 
schemes (RS) within the dimensional- regularization prescription (see, e.g., Refs. JTT|, |T2"|] ) . 

The two variants of dimensional regularization that are mostly used in one-loop com- 
putations are conventional dimensional regularization (CDR) and dimensional reduction 
(DR). In both schemes one considers 2 helicity states for spin-| Dirac fermions. The es- 
sential difference between the two schemes regards the number of the helicity states of the 
gluons in the loop. The gluon has d — 2 helicities in CDR and 2 helicities in DR. Since 
the number of helicity states of gluons and Dirac fermions is the same in DR, this scheme 
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preserves supersymmetric Ward identities. Within each scheme one can still choose the ex- 
ternal particles (their momenta and helicities) in the amplitudes to be either <i-dimensional 
or 4-dimensional. At one-loop order, these choices lead to differences of 0(e), which do not 
have any effect on the results presented in this paper. 

Note that ultraviolet and infrared divergences behave differently with respect to di- 
mensional-regularization prescriptions. The ultraviolet RS dependence can ultimately be 
removed by a proper redefinition of the renormalized running coupling. The infrared RS 
dependence instead leads to contributions (which are not vanishing for e — > 0) that depend 
on the specific amplitude and that cannot be reabsorbed by an overall (i.e., independent 
of the amplitude) redefinition of the renormalized coupling. These features, which were 
explicitly pointed out in Ref . jl]J , will be discussed in detail in the following sections. Note 
also that, in the calculation of physical quantities, the RS dependence of loop amplitudes 
has to be consistently matched to that of tree amplitudes. This issue is discussed on quite 



a general basis in Ref. [12]. 



2.2 Partially renormalized amplitudes 



We consider amplitudes A m that involve m external coloured particles (gluons, massless 
and massive quarks, gluinos and squarks) with momenta pi, . . . , p m , masses mi, ... , m m and 
an arbitrary number and type of colourless particles (photons, leptons, vector bosons, etc.). 
We always consider the amplitudes in the crossing-symmetric, but unphysical channel when 
all particle momenta are outgoing. The amplitudes are denoted by A m (pi, mi, • • • ,p m , m m ) 
(or, shortly, A m ({Pi, mj)), and the dependence on the momenta and quantum numbers of 
the colourless particles is always understood. 

For amplitudes of processes involving massive particles the MS subtraction scheme is 
not always used to perform charge (coupling) renormalization (see, e.g., Ref. [0). To leave 
the freedom of choosing a favourite charge-renormalization scheme, we find it convenient to 
use mass-renormalized, but charge-unrenormalized amplitudes. Thus, in the perturbative 
expansion 



A m (g s ,fi 2 ; {pi, mi}) 



9sH 

47T 



g s stands for the bare strong coupling and m, are the renormalized mass parameters (i.e., 
masses and related parameters such as those which appear in Yukawa couplings). The 
renormalized masses are obtained from the bare masses mf 1 ' by the replacement mf' — > 
rrii = mf + g B Srrii, so that the (ultraviolet-divergent) mass renormalization constants 
Srrii are implicitly contained in A$ ■ In Eq. (|l|) q is a non- negative integer, and /x is the 
dimensional- regularization scale. 

Equation ([TJ) fixes the normalization of the tree- level, A$, and one-loop, A$, coefficient 
amplitudes that we use in the rest of the paper[| Although it is not explicitly denoted in 

^Precisely speaking, Am is not necessarily a tree amplitude, but rather the lowest-order amplitude for 
a given process; A m is the corresponding NLO correction. For instance, in the case of gg — * 77, AW 
involves a quark loop. 
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Eq. (HD, A^ and A^ are both dependent on the RS. 



2.3 Colour space 

We shall present the singular structure of the QCD and SUSY QCD amplitudes directly in 
colour space. In particular, we use the same notation as in Ref. ||. 

The colour indices of the m partons in the amplitude A m are generically denoted by 
Ci, . . . , c m : Cj = {a} = 1, . . . , N% — 1 for particles in the adjoint representation (gluons, 
gluinos) and q = {a} = 1, . . . , N c for particles in the fundamental representation (quarks, 
squarks and their antiparticles) of the gauge group. We formally introduce an orthogonal 
basis of unit vectors {\ci, . . . ,c m )} in the m-parton colour space, in such a way that the 
colour amplitude can be written as follows: 

A c £- ,Cm {pi,m u . . . ,p m ,m m ) = (ci, . . . ,Cm | Amipum,!, . . .,p m ,m m )) . (2) 

Thus \A m (pi,mi, . . . ,p m ,m m )) is an abstract vector in colour space, and the square am- 
plitude summed over colours is 

\A m ({pi,mi})\ 2 = (A m ({Pi,mi}) \ Am({pi, m»})) . (3) 



Colour interactions at the QCD vertices are represented by associating a colour charge 
Tj with the emission of a gluon from each parton i. The colour charge Tj = {T"} is a 
vector with respect to the colour indices a of the emitted gluon and an SU(N C ) matrix 
with respect to the colour indices of the parton i. More precisely, its action onto the colour 
space is defined by 

(ci, . . . , q, . . . , c m I If | 6i, . . . , hi, . . . , b m ) = 6 cibl . . . T c ^ . . . S Cmbm , (4) 

where T^ b is the colour-charge matrix in the representation of the final-state emitting par- 
ticle i, i.e. T°jj = if cab if z is a gluon or a gluino, T®p = t a afi if i is a (s)quark, and T£p = —tp a 
if i is an anti(s)quark. The colour-charge operator of an initial-state parton is defined by 
crossing symmetry, that is by (Tj)^ = — t% a if % is an initial-state (s)quark and (Tj)^ = t a a ^ 
if i is an initial-state anti(s)quark. 

In this notation, each vector \Am(pi,mi, . . . ,p m ,m m )) is a colour singlet, so colour 
conservation is simply 

m 

^T,|^l m ) = 0. (5) 

i=l 

The colour-charge algebra for the product (Tj) a (Tj) a = Tj ■ Tj is: 

Ti ■ Tj — Tj • Ti if i^j; Tl = Ci, (6) 

where Ci is the quadratic Casimir operator in the representation of particle i, and we have 
C F = T R {Nl - l)j N c = (N 2 C - 1)/(2AQ in the fundamental and C A — 2 T R N C = N c in the 
adjoint representation, i.e. we are using the customary normalization Tr = 1/2. 
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We remind the reader that, in the cases of amplitudes with m = 2 and m = 3 coloured 
partons, the colour-charge algebra can always be recast in a fully factorized form in terms 
of Casimir operators of the m partons (see the Appendix A of the second paper in Ref. ||). 

Note that in the following sections we always refer to one-loop QCD corrections to 
scattering amplitudes. Nonetheless, most of the results can straighforwardly be used for 
the case of one-loop QED corrections. To this purpose it is sufficient to replace the colour 
couplings g a Ti by the electric couplings gei, where g is the electric charge of the electron 
and Ci is the charge of the parton i in units of the electron charge. In terms of colour 
factors, this implies the replacements Cp — > ef, T R — > 1 and Ca — > 0. 

3 Singular behaviour at one-loop order with massive 
partons 

3.1 QCD amplitudes 

In this subsection we present our results on the singular behaviour of QCD amplitudes 
at one-loop order. In the massless case the one-loop coefficient subamplitude A$ has 
double and single poles in e that can be obtained by a process-independent factorization 
formula [§-0. Similar poles, although with different coefficients, appear if the amplitude 
involves massive partons. We find that these singularities are still universal, so that the 
factorization formula for the massless amplitudes can be generalized to the massive case. 

The general factorization formula for the one-loop coefficient subamplitude A$ isffl 

(7) 

All the e-poles are included in the factor I, so the remaining contributions on the right- 
hand side can be safely expanded in e for e — > 0. Moreover, the factor I includes also the 
constant (when e — > 0) terms related to the RS dependence. Thus, the contribution A$ &n 
is not only finite, but it is also RS-independent. 

These features are shared by massless and massive amplitudes. In the case of massive 
quarks, moreover, our factorization formula has additional important properties related 
to logarithmically-enhanced contributions. In the limit of one or more vanishing masses 
mj, the subamplitude A$([J?; {pi, wii}) contains logarithmic terms of the type ln 2 mj and 
In rrti that become singular. It also contains constant terms that originate from the non- 
commutativity of the limits rrii — > and e — * 0. We are able to embody all these logarithmic 
and constant terms[J] in the factor i". Thus, in the limit — ► the finite contribution 
(• • • iVii m «; • • •) tends smoothly to the finite contribution A$ &n (. . . ,Pi,rrii = 0, . . .) 

^Here and in the following, the labels R.S. explicitly denote the RS dependence of the various quantities. 

"in the small-mass limit, the mass m t replaces dimensional regularization as regulator of collinear 
singularities. In this sense, we can say that we can control the ensuing regularization-scheme dependence 
including finite (when wii — > 0) terms. 
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of the corresponding amplitude in the theory where the parton i is massless: 

lim ^ )fin (ju 2 ;pi,mi, . . . } p i} m i} . . .) = A^ Rn (/j, 2 ; pi, m u . . .,jjj,m< = 0, . . .) . (8) 

Note that Eq. (H) is valid independent of the actual definition of the renormalized mass 
rrii (or of the related Yukawa couplings). We can use either the pole- mass definition or the 
MS definition, because the terms In rrii originating from the different definitions are always 
suppressed by the corresponding mass factor as rriilnrrii. 

The factorization formula (|7|) and the property in Eq. (H) can be used to check the calcu- 
lation of the massive amplitude by comparing it to the corresponding massless calculation. 
In the asymptotic regime where rrii is much smaller than any of the relevant kinematic 
invariants Q, these equations can also be used to directly obtain (apart from corrections 
of Oirrii/Q)) the one- loop massive amplitude from the corresponding massless amplitude, 
without explicitly computing the former. 

The first term on the right-hand side of Eq. (0) has a factorized structure in colour 
space. The singular dependence (poles in e and logarithms in rrii) is embodied in the factor 
J^f that acts as a colour-charge operator onto the colour vector |.A$P) RS . Note that both 
factors are RS-dependent. In particular, the product of the RS-dependent terms of 0(e) in 
A$} and double poles 1/e 2 in J produces, in general, an RS dependence of A$ that begins 
at 0(l/e). 

The explicit expression for I m in terms of the colour charges of the m partons is the 
following: 

r(l-e)\ q 2 \ e 



/i 2 ; {p z , m,}) = q - ™ - ft 



+ E T rT k 



v£W;^-^) + - (- 



Vjk \ e 2 




3=1 



The first term on the right-hand side of Eq. (||) contains the ultraviolet divergences, to 
be removed by the renormalization of the bare coupling g s . It is proportional to q, which is 
the overall power of g s in Eq. ([]]), and (3q is the first coefficient of the QCD beta function 

fo = ^C A -^T R (N f + N F ), (10) 

where Nf and Np are the numbers of massless and massive quark flavours, respectively. 
The constant coefficient /3q a parametrizes the ultraviolet RS dependence. Setting /3q DR = 
by definition in the CDR scheme, its corresponding value in the DR scheme is |nj (see also 
Ref. 0) 

$ K = \C A . (11) 

The second and third terms on the right-hand side of Eq. @ have an infrared origin. 
We have determined them by exploiting the fact that, in NLO calculations of infrared- 
safe cross sections, the infrared singularities of the one-loop amplitudes are cancelled by 
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their counterpart in the real-emission contribution. The latter has been computed by 
extending the dipole subtraction formalism || |9| to the case of massive partons. The 
infrared contribution that leads to colour correlations proportional to TV,- ■ is produced 
by soft and collinear singularities. The infrared contributions r^ a are produced by either 
collinear (but not soft) or soft (but not collinear) singularities. 

The singular function Vj^ that controls colour correlations is symmetric with respect 
to j <-> k and depends on the Lorentz invariant Sjk = 2pj ■ Pk and on the parton masses. In 
particular, it depends on the the relative velocity Vjk of particles j and k: 



ijfn 2 
~\2 



V jk 

Its explicit expression for non-vanishing masses rrij and nik is 



,1-—^. (12) 



V<?>(« m,, ,n k ; e) = I ±- Ini^tt - i(ln^ + - £ , (13) 

2e v jk 1 + Vjk 4 y 



while for one or two vanishing masses we find 



V«ffe;0,0; e ) = i. (15) 



The functions depend on the flavour of the parton j and on the parton masses. In 
the case of gluons and massless quarks (antiquarks) we have 

1 2 Nf m 2 

r-Gu,m w ;e) = - lg - 7™ - - T R £ In -f , 
e o F=1 /i 

rr(^0;e) = \l q -%\ 
while for massive quarks (antiquarks) we find 

CI 7T? \ TV? 

e 11 j 11 

The flavour coefficients jj in Eqs. (]I6|) and flTTD are 

7i^,j=§C,, 7g = H CA _^ Ti?iV/ . (19) 

The coefficients 7j a parametrize the finite (for e — > 0) contributions related to the RS de- 
pendence of the one-loop amplitudes with external massless partons (the massive-quark 
function in Eq. (plf ) does not depend on the actual version of dimensional regularization 
used in the calculation). The transition coefficients 7j s that relate the RS mostly used in 
one-loop computations were first calculated in Ref. [[11]. They are given by 



(16) 
(17) 



1 1, m 

e 2 /i z 



Note that the dependence on the RS of the ultraviolet and infrared contributions is 
different. The ultraviolet contribution (3^ & in Eq. (|9|) is proportional to the overall power 
q of g s that controls the amplitude (see Eq. (jH)) and it can ultimately be reabsorbed by a 
process-independent redefinition of the renormalized coupling. The infrared contributions 
7j" a to Eq. (|D instead depend on the number and flavour of the external massless partons in 
the amplitude. In the calculation of physical (infrared-finite) quantities this dependence has 
to be cancelled by computing the corresponding real-emission contributions in a consistent 
manner, that is, by using the same dimensional-regularization prescription as in one-loop 
amplitudes |12| . 



Some comments on the structure of these results, in particular about the e poles and 
the mass logarithms, are appropriate. 



If there are only massless partons, our result agrees with those in Refs. |J-|6[. The 
structure of the e poles for massive quarks agrees with the results of Ref. || and with 
the QED case considered in Ref. |§. 

The double poles 1/e 2 in Eq. (|7|) are factorized completely and not only in colour 
space. More precisely, there is a contribution — T 2 /e 2 from each external massless 
partons j. The simplest way to see that is to expand Eq. (^) in powers of e and then 
use the colour conservation relation (|), i.e. T fc = — Tj. One obtains the result 

I m (e, /i 2 ; {pi, m,}) = \ll T r T k + 0(l/e) = -\ £ T 2 + 0(l/e) (21) 

j € k^j € J 

mj =0 m,j — 

that explicitly shows the absence of colour correlations at 0(l/e 2 ). Nonetheless, single 
poles 1/e are both colour- and velocity-correlated. 



The term proportional to 1/e in Eq. ([13]) is familiar from QED bremsstrahlung, and 
the imaginary part in Eq. (|9|) is the corresponding Coulomb phase. 



Since we have 

,2,^,2 

(22) 



1 ~ Vjk rntmj. 

l + v jk s) 



5 ifc 

Eqs. (|I3"D, (P|) and (|i5|) are related by the following formal correspondence 

1 , m 2 1 , „ m 2 7T 2 m ( m 2 , \ 1 , . 

— ln^--ln 2 ^ + O -r—2-r < >— (23) 

2e \s jk \ 4 \s jk \ 12 \\Sjk\J 2e 2 

between mass logarithms in the massless limit and 1/e 2 poles in the massless theory. 

• The various constant (for e — > 0) terms in Eqs. (|i"3|)-([l8|) are important to guarantee 
the smooth limit in Eq. (|8|). One can always include additional finite terms in J, 
provided they are smooth in the massless limit. We chose to include the terms that 
are proportional to 7i 2 Q(sjk) in Eq. (Q). 

We also add some other comments on the origin of the last two terms on the right-hand 
side of Eq. (§). The contributions r^ a can be related to the Altarelli-Parisi splitting 
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functions (see Sect. |J). As already mentioned, we have evaluated the colour-correlation 
term in Eq. ([)]) by computing the corresponding bremsstrahlung contribution. This directly 
gives the real part of this term, namely the function vj^ . To obtain the corresponding 
imaginary part in the square bracket on the right-hand side of Eq. (|j), we have exploited 
the fact that the singular (real and imaginary) part of the colour-correlation term is directly 
proportional to the the following three-point function: 

T T [ d d a % PjPk 

J k J q q 2 + iO [( Pj + q) 2 - m 2 + iO] [(p k - q) 2 - m 2 + iO] ' K ' 

The contribution of the gluon pole, 

-^—:^2n5 + (q 2 ), (25) 
q z + «(J 

gives Vj k , while the contribution from the poles in the massive propagators, 

— 1 2— — - -> -47r 2 5+((p j + qf - m 2 ) 8+((p h - q) 2 - m 2 ) , 

(26) 



gives the corresponding imaginary part (see, e.g., Ref. fig]). Note that in the massive case 
the imaginary part is more involved than in the massless case, where it is obtained from the 
real part by the simple analytic continuation \n(sj k ) — > ln(— Sj k — iO) = In \sjk\ — iTrQ(sjk) 
of its overall factor (fi 2 /sjkY i n Eq. (0). 

According to our definition of the insertion operator i" in Eq. (0), the finite one-loop 
contribution A$ &n still depends on the dimensional-regularization scale \x. This depen- 
dence is nonetheless simple, because it is embodied in a contribution proportional to In \i. 
The coefficient of this single-logarithmic contribution is given by 



d 



2 



/^^ )fi V; OW) = U f - E [l 3 - ^T R N F 5 jg \j A^({ Pl ,m t }) , (27) 

Note that /i should not be confused with the renormalization scale. In particular, if renor- 
malized masses (and related Yukawa couplings) do not correspond to the pole-mass defi- 
nition, but they are MS-scheme running masses at the renormalization scale /xr, the finite 
contribution *4^ fin can explicitly contain terms of the type mjln(mi/ /xr). 



3.2 SUSY QCD amplitudes 



All the results on QCD amplitudes in Sect. |3J] can be extended to SUSY QCD processes 



A preliminary discussion on the RS and on the small-mass limit is nonetheless necessary. 

For a supersymmetric theory CDR (in contrast to DR) is not a consistent RS, because 
the mismatch between the d — 2 = 2(1 — e) transverse degrees of freedom of the gluons and 
the 2 transverse degrees of freedom of the gluinos violates supersymmetric Ward identities. 
In particular, in the case of on-shell one-loop amplitudes, this leads to violation of the tree- 
level identity g s = g s , between the gluon (gauge) coupling g s and the qqg- Yukawa coupling 
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g s . Nevertheless, CDR can also be used in SUSY calculations, because supersymmetry can 
be restored by introducing a proper counterterm ||T6|,|T7||. More precisely, we can still use 



the notation of Eq. (|TJ) (where only g s appears) in any RS, provided we implement the 
following RS-dependent relation between the two couplings at one-loop order: 



9s \ 2 »M 



(28) 



where 



7 =v, 7 =^C A --C F . (29) 



As discussed in detail in Sect. |3.1| , our factorization formula for one-loop QCD ampli- 
tudes smoothly interpolates between the cases of massless and massive quarks. Analogous 
results, with smooth behaviour with respect to sparticle masses, can be obtained for SUSY 
QCD. However, in this section we are not going to present such general results because 
of the following reasons. On one hand, the case of exactly massless SUSY partons has no 
practical interest. One the other hand, owing to the proliferation of masses in the SUSY 
particle spectrum, the control of logarithmically-enhanced terms produced by the smallness 
of the mass of some SUSY partons would require a quite involved presentation and detailed 
discussions of the various possible cases related to different scenarios of mass hierarchies. 

For the sake of simplicity, we thus limit ourselves to considering the case of finite (non- 
vanishing) squark and gluino masses, rrig and m~ g . The factorization formula given below 
includes all the divergent e-poles that appear in this case, and some related logarithmic 
terms. 

The one-loop amplitudes for SUSY QCD processes have the same singularity structure 
as given by Eqs. (|7]) and @. The coefficient /3 in Eq. (|||) should simply be replaced by 
/^ USY , the first coefficient of the SUSY QCD beta function: 

^susy = H ^ _ 4 Tr{ ^ + Np) _2 TrNs _ 2 Ca ^ (3Q) 

where N s is the number of chiral squark pairs q = (g#, q F ) (N s — Nf + N F in the fully 
supersymmetric theory). Moreover, the sums over the flavour indices j, k in Eq. ([|) run over 
all parton species (j, k = q,q,g,q,g). Correspondingly, we have to introduce the flavour 
functions Tj of the external massive gluinos and squarks: 



T j (iJ,m r ,e)=T 2 J - - hi ^ - 2 \ +lj hy-£, J = ~g,q, (31) 

where 

H = \Ca, % = 2C f . (32) 



As discussed above, the extension of Eq. @ to SUSY QCD by using Eqs. (|3~0D-(j32l) 
does not include in the insertion operator i" all the terms that are logarithmically-enhanced 
in the limit of small squark and gluino masses. Note however, that the finite contribu- 
tion A^^ n (fj, 2 ; {pi^rrii}) of the one-loop SUSY QCD amplitude still fulfils the property 
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in Eq. (||D, provided the massless limit m ; 
(i = q,q). 



is restricted to quarks and/or antiquarks 



The /x-dependence of A$ An (fJ?; {pi, rrii}) is still given by Eq. (|27f) , with the replacement 
Po —> Pq VSY an d the inclusions of the flavour coefficients Jj=g t q in Eq. fl32|). 



Relation of the singular terms to the Altarelli-Parisi 
splitting functions 



In the case of massless QCD, it is known that the flavour functions r^ a of Eq. (|S|) 
are related to the Altarelli-Parisi splitting functions. In this section we first recall this 
correspondence and then we sketch how it can be extended to the massive case. 

Let us consider, for instance, the collinear splitting of a massless quark. The correspond- 
ing Altarelli-Parisi splitting function P^ q & (z; e) in the dimensionally-regularized theory can 
be written as 



3RS / 
qq \ 



1 + z 2 



+ 



hi 



2C f 



(33) 



where h^ & (e) is the number of gluon helicity states in the corresponding RS (see Eq. ([42])). 
On the right-hand side we have isolated the contribution that becomes singular in the soft 
region z — > 1 from the remaining regular part P g ^ cg ) RS . 

The soft contribution^ has to be combined with the soft terms produced by non- 
collinear (large- angle) gluon radiation (see Eq. (|2~4])). The terms that are singular both in 
the soft and collinear regimes lead to the double-pole singularity of the colour-correlation 
function in Eq. (|i~5"l). In the massive case kinematics is different from the massless case, 



and the contribution of these soft terms is embodied in the functions vfk\ s jk'i mj, n^k] e), 



which control the entire colour-correlation part of Eq. (||]). 



The term p^ e s) Ra can only lead to collinear singularities (single poles and constants) 
and it is related to the function Tg & in Eq. (|TTD, including its RS dependence 



7? 



~CDR 



0) 



dz - 
o e 



p(reg)F 



qq 



p(reg)CDR, 



(34) 
(35) 



Analogous equations relate the other flavour coefficients jj, 7j*- a in Eqs. (|T9|), (|20|) , (|32|) to 
the corresponding terms P^ eg ' )Ra of the Altarelli-Parisi splitting functions. 

The relation between the functions r RS and the splitting functions can be extended from 
the massless to the massive case, provided we properly take into account the corresponding 
dynamics and kinematics differences. 



**Analogous soft terms, which become singular when z — > 0, appear in the Altarelli-Parisi functions of 
other splitting processes (see the equations in the Appendix). 
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The dynamics of the splitting processes of massive partons can be described by gener- 
alizing the collinear limit to the quasi- collinear limit. Let us consider a generic tree-level 
amplitude A.m+i({Pi, m^}) with m + 1 external partons. The limit when an internal parent 
parton (labelled by (jk)) decays quasi-collinearly in two external partons j and k is defined 
by 

P? - , Vl - (1 - z)p>* , P 2 = m 2 m , (36) 



with the constraint 



mj ml m 2 w 



Pj ■ Pk, m~, m&, rrifjk) ^ at fixed ratios — , , — - — - . (37) 

Pj ■ Pk Pj ■ Pk Pj ■ Pk 

The quasi-collinear limit obviously differs from the collinear limit because the splitting 
partons are massive. However, the key difference between the two limits is given by the 
constraint that the on-shell masses squared have to be kept of the same order as the 
invariant mass (pj + Pk) 2 , while the latter become small. 

It can be shown that in the quasi-collinear limit the tree-level squared amplitude fulfils 
a factorization formula similar to the analogous formula for the collinear limit. We have 

l-^il 2 - l^l 2 fp + T)f s m 2 p?AM * to 2 }) . (38) 

{Pj T Pk) m (jk) 

where the m-parton matrix element on the right-hand side is obtained from ^4^+i by 
replacing the partons j and k with the single parent parton of momentum p. The function 
P(jk),j( z: i e i {/•*?}) generalizes the customary <i-dimensional Altarelli-Parisi splitting function 
to the case of quasi-collinear splitting. It has the usual dependence on e and on the 
longitutinal-momentum fraction z, plus an additional dependence on the parton masses 
indicated by the variables 

2 

[Pj + Pk) 2 ~ m'jk) 

where (in general) mf stands for any quadratic combination of the masses of the partons 
involved in the splitting process (mf = m 2 , m 2 ., m 2 j k ^,mjmk, ■ ■ ■ and so forth). 

For instance, in the case of the quasi-collinear splitting q — > q + g of a massive quark 
(p 2 = m 2 , p 2 = 0), the analogue of the massless splitting function in Eq. (|33|) is (see 
Eq. (H) and Refs. @,1): 



^ [ Z;e; 2p q pJ- l-z + 



2 

777 
PqPg\ 



(40) 



The 'regular' part of the splitting function is now given by the contribution in the square 
bracket. It contains a term that explicitly depends on m q , in addition to the regular part 
that appears in the massless case (see Eq. (pB])). 



Comparing the regular parts of the splitting functions in Eqs. (|33|) and (fEj), we can 
understand the difference between the the massless and massive functions T q in Eqs. (pp 
and fll8|) . This difference amounts to the correspondence 

'"" A ■ (41) 
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In both the massless and massive cases, 7 9 is given by Eq. fl34"|). The coefficient in front of 
7 g is obtained by performing the integration of the propagator factor l/p q p g (see Eq. fl3"8|)) 
over the relative angle 9 qg between the quark and the gluon. In the massless case this 
integration is singular when qg — > and it leads to the single pole 1/e on the left-hand 
side. In the massive case, the integration is kinematically regularized at a cutoff value 
9 qg ^ m q j E q and it leads to the logarithmic behaviour of the first term on the right-hand 
side. The contribution inside the round bracket on the right-hand side is instead produced 
by the term of Eq. (|40|) that explicitly depends on m q . Here, the angular integration, which 
is dominated by the region 6 qg ^m q /E q , produces a constant term, while the integration 
over the gluon energy E g ~ (1 — z)E q is divergent in the soft region z — > 1 and it leads to 
the single pole 1/e in front of the round bracket. 

The generalized Altarelli-Parisi functions for the quasi-collinear limit of the other 
splitting processes in QCD and SUSY QCD are listed in the Appendix. The expressions 
of the flavour functions in Eq. (16) and Eq. ( j3l"l) are related to the splitting functions in 



Eqs. (M), (H) and Eqs. (M), (EHT), respectively. 



5 Summary 

In this paper we have discussed the singular behaviour of on-shell QCD and SUSY QCD 
amplitudes at one-loop order in the presence of massive particles. The complete structure 
of the ultraviolet and infrared singularities is described by the colour-space factorization 
formula given in Sect. ^. The factorization formula is universal, i.e. valid for any amplitude, 
and the singular factors only depend on the flavours and momenta of the coloured external 
legs. Moreover, the factorization formula is given in such a form that the corresponding 
formula for massless QCD partons is smoothly recovered by simply letting the masses 
approach to zero. 

Our factorization formula can be useful both to check explicit evaluations of one-loop 
amplitudes and to organize their calculations in terms of divergent parts and finite remain- 
ders. Furthermore, in the asymptotic regime where the parton masses are much smaller 
than any of the relevant kinematic invariants, the formula can also be used to directly 
obtain (apart from vanishing corrections when the masses tend to zero) the one-loop mas- 
sive amplitude from the corresponding massless amplitude, without explicitly computing 
the former. In the general context of NLO calculations of jet observables, our one-loop 
results are useful for setting up the integration of tree-level amplitudes in such a way as to 
construct process-independent techniques for infrared cancellations. 
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Appendix: Quasi-collinear dynamics and generalized 
Altarelli— Parisi splitting functions 



In this appendix we list Altarelli-Parisi splitting functions in QCD and SUSY QCD, supple- 
mented by the mass terms that are relevant in the quasi-collinear limit discussed in Sect. f| 
(see Eqs. (BB)-(BO)). In the expressions given below all the mass terms can be parametrized 



by the variable fj, 2 k 



ml 



+ ml)/[(pj+p k y - 



m 



Cj'*)J 



We include the dependence on the 



RS, which is parametrized by the corresponding number h^ & (e) of gluon polarizations, 



hf DK = d-2 



2e 



h 



DR 



(42) 



Note that our expressions refer to Altarelli-Parisi splitting functions after average over the 
azimuthal angle identified by the collinear direction. The azimuthally averaged functions 
are relevant to discuss the RS dependence [jT2| and, in general, they do not coincide with 
the splitting functions averaged over the polarizations of the parent parton. Since we 
are considering the splitting process (jk) —>■ j k, the splitting functions fulfil the obvious 
symmetry relation Pffi )k {z; e; fi 2 kj ) = P$ s k)j {l - z; e; fi 2 jk ). 

Quarks and gluons: 



Pqq i Z i £■> f^qq) -^~J 



R 



2C, 



'-^ z + \h?{l-z)-2»l 
1 " (2*(1 -z)- Mqq 



d-2 

z 1 — z 
+ + 



h °*-z(l-z) 



d-2 



(43) 
(44) 
(45) 



Gluinos and gluons: 



Pfi( z \ e ! ^gg) 

P K ~ & (z-e- u 2 ~ 

1 gg \"i c i rgg) 



C A 
C a 



99 



1 - (2z(l -Z)-LL 

d - 9. V v > p 



(46) 
(47) 



Squarks and gluons: 



gq V*) c ) H'qqJ - n j _ ^ 



c F 

Tr 



2z 



2(JL 



2z{\ - z) 



(48) 
(49) 



Equations 
quarks. 



i-(|4*9|) are valid for the superpartners of both the left- and the right-chirality 



Neglecting the mass terms, which are proportional to fi 2 k , the splitting functions coincide 
with those reported in Ref. [17[] for the CDR scheme. We do not consider the splitting 
functions produced by the Yukawa coupling qgq, because, as long as the gluino and squark 
masses are finite, they do not produce any singular terms when e — > 0. 
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